Markovian feedback to control continuous variable entanglement 
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We present a model to realize quantum feedback control of continuous variable entanglement. It 
consists of two interacting bosonic modes subject to amplitude damping and achieving entangied 
Gaussian steady state. The possibility to greatly improve the degree of entanglement by means of 
Markovian (direct) feedback is then shown. 
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Entanglement is one of the most puzzling features of 
quantum mechanics that has given rise to long debatés 
on foundational aspects since the seminal papers on the 
subject 0,0- Recently, it has been also recognized as 
a valuable resource for quantum information processing 
0- Thus, its generation and control has became a pri- 
mary task to accomplish. So far much attention has been 
devoted to generation and characterization of entangle- 
ment, while for its control one usually refer s to error cor- 
rection or distillation procedures Q ■ 

Feedback provides an intuitive way to control a system 
state, hence entanglement since it is a system state pecu- 
liarity. A theory of quantum-limited feedback has been 
introduced by Wiseman and Milburn jj] leading to rel- 
evant experimental achievements yj- It is based on the 
direct (immediate) use of the measurement results to al- 
ter the system state, hence the name Markovian feedback. 
It turns out to be particularly useful in realizing squeezed 
states [f|. In fact, monitoring a system's observable will 
conditionally squeeze its variance. Unconditional squeez- 
ing is then obtained by using the measurement results to 
continuously drive the system into the desired, determin- 
istic, squeezed state. If such arguments are applied to a 
multiparties system in case of nonlocal measurements, 
the realized squeezing is related to entanglement among 
the parties Q. 

However, the crucial point is to see whether local mea- 
surements followed by feedback action sufïice to control 
entanglement (because local measurements condition the 
system to separable states). We shall show that it is in- 
deed possible. By referring to Fig^ we can consider two 
interacting subsystems SI and 52 each one loosing infor- 
mation on its own environment El and El. Then, mon- 
itoring such environments separately will give outcomes' 
currents that can be joined and exploited for feedback 
action. 

Recent results on finite dimensional systems show that 
such feedback scheme is helpful in enforcing entangle- 
ment However, the original gedanken experiment 
proposal of Ref.0 refers to infmite dimensional systems; 
moreover continuous variables offer a lot of advantages 
in information processing @. We henceforth focus our 
attention on the possibility to control the degree of en- 
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FIG. 1: The studied scheme. Two subsystems S'l and 5*2 
interact (thik dashed line). Each one is coherently driven 
(left-right arrow and right-left arrow) and loses information 
in its own environment, El and E2. The environments are 
separately monitored and then a joint current is exploited for 
a feedback action on the driving fields (dashed lines). 



tanglement for a two-party quantum Gaussian state. A 
paradigm of the model of Ref . 0] is represented by a non- 
degenerate parametric oscillator wherc two bosonic 
modes a\ and a 2 interact through a quadratic Hamilto- 
nian of the type (h = 1) 
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with x £ K the coupling strength. In addition we con- 
sider the two modes driven by a Hamiltonian 



H 



drive 



(^a\ + <4) + a * ( a i + a 2) 



(2) 



where a£Cis the driving amplitude (assumed equal for 
both modes). Hence, the total Hamiltonian results 
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Now, suppose that the system looses information on 
each mode through an amplitude damping channel at 
rate k (hereafter we consider k = 1 so that all physical 
quantities become dimensionless). Then, by introducing 
the modes quadratures 
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the system dynamics can easily be described by the quan- 
tum Langevin equations 
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where Xj and 3^- (j = 1, 2) are hermitian operators rep- 
resenting the vacuum noise entering the system. They 
have equal time correlations 



(XjXf.) 
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By monitoring each mode environment it would be pos- 
sible to measure a quadrature for each subsystem, e.g. 
by using homodyne detection 01 • Suppose this is done 
with an overall efficiency ry S [0, 1] (i.e., 77 accounts for 
the detectors efSciency, the fraction of the field b eing 
measured, etc), then we can write local currents like |l·l| 



Ij(t) = vXj + VvWj, 3 = M : 



(11) 



where Wj (j = 1, 2) are hermitian operators representing 
the vacuum noise affecting the currents. They have equal 
time correlations 
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We can now combino the currents Ij to get the following 
joint current 



I(t) =r l {X 1 - Xa) +^j(Wi-W 2 ) 



(13) 



through which it would be possible to gain information 
about the quantity X\ — X 2 . Its variance in absence of 
environmental effects and for increasing value of x, tends 
to zero |l0) . corresponding to the maximally entangled 
state of the type discussed in Ref.0- Thus, the fact that 
it goes below 1 (vacuum fluctuations) can be roughly con- 
sidered as signature of entanglement. In such a sense, we 
may assume the current in Eq. (|13f) giving us information 
about the system state entanglement. 

The steady state solution of Eqs.(|5j)-|0) results 



(Xi) ss = (X 2 ) ss = -iV2 
(Yi)„ = (Y 2 ) ss =-V2j 
which is stable for x < 1/2. 
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Subtracting Eq.JTJ) from J5J we get a quantum 
Lan gevi n equation resembling a Ornstein-Uhlenbeck pro- 
cess [13 ■ It can be easily solved to get the steady state 
variance 



((X 1 -X 2 ) 2 )-((X 1 ) SS -(X 2 ) SS ) 5 
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l + 2x 



(16) 



We see that the variance Ijl6(l goes below 1 (vacuum fluc- 
tuations) as soon as \ > 0, thus we infer the presence of 
entanglement in the steady state. However, since it must 
be x < 1/2, the variance (|l(jf) is limited from below by 
1/2. By converse, the amount of entanglement will be 
limited. That is, the damping channel degrades the sys- 
tem state preventing it to become maximally entangled 
like that of Ref.0. 

We can now think to control such a process, hence to 
control entanglement, by using feedback action accord- 
ingly to the information gained about X\ — X 2 . To this 
end we consider an additional Hamiltonian proportional 
to the current l|13l) and driving the operator Y\ — Y 2 , 



H fb = -I(t-r) (Y 1 -Y 2 ) 
V 



(17) 



Here A represents the feedback strength and r the feed- 
back loop delay time. The above choice of the driving is 
motivated by the fact that we want to affect the variance 
of Xi — X 2 . The feedback action 117|) can be realized 
by a modulation of the driving fields accordingly to the 
current /. 

Adding the Hamiltonian 1|17|) to the Langevin equa- 
tions ©-© means to add, for a genèric operator O, the 
tcrm 



O fb = - / drG(T)I(t - r)A [Y, -Y 2 ,0] , (18) 
V Jo 

where G is the feedback response function Q. Practi- 
cally, it accounts for the feedback loop delay time, how- 
ever hereafter we assume it negligibly small, hence G can 
be considered as a Dirac delta function. 

When inserting Eq. i(T%|) into Eqs.JSJ)-©, one has to 
also account for the noise carried inside the system by 
the current l|13|) . In presence of imperfect detection the 
current may be non trivially related with the input noise. 



put : 
.14] 



The following correlations hold at equal times 

ÇWjXk) = (x j W k ) = ^S ]k , (19) 

(Wjy k ) = -(y,w k ) = i^-5 Jk , j, k = 1,2 . (20) 

It is immediate to see that in the perfect efficiency case 
(77 = 1) one can identify the feedback noise with the input 
noise, while in the opposite case (77 = 0) the two noises 
are uncorrelated (as it can be easily expected since in 
such a case the feedback noise has nothing to do with 
the vacuum input noise). 
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It is worth noting that the feedback action l|18|l does 
not affect the steady state vàlues (|14|) - (|15ll . Then we can 
rewrite Eqs. (EJ- (ISJ i n presence of feedback, for only the 
steady state fluctuations 



xj — Xj (Xj) ss , 
Vj = Yj-(Yj) ss , i = 1.2. 
Let us introduce the operator vectors 
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then, the feedback modified Langevin equations read, in 
a compact way, 



where 
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The symmetric noise correlation matrix is 



N 



1 



± ((nn T ) + (nn T ) T ) 



° , * 

-A — — 

n 





-A — — o\ 
v ' 



i + A+^ 



(26) 



therefore the stationary symmetric correlation matrix 



1 



t ((vv T ) + (vv T f) , 



(27) 



2 

can be derived from the typical relation for Ornstein- 
Uhlcnbeck-like processes |l3j 

MT + TM T = N . (28) 



For the following purposes, we write T in terms of its 
2x2 submatrices 



7 a 

T 

a 7 



where the matrix elements, by virtue of Eq, 

(1/2) - a + (i - 2 X ) (x + xy v ) 
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1 - 4A + 8A X - 4 X 2 
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2 Vl-4x 2 
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and 
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X-A-(1-2 X ) (A + AVt?) 

1 - 4A + 8A X - 4x 2 
0-21=0, 
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Since the steady state is Gaussian it is completely char- 
acterized by the correlation matrix 1)2 7[l. Then, its degree 
of entanglement can be quantyified by means of the log- 
arithmic negativity jlól| 



L : 



log(2C) 




ifC < i 

otherwise 



(36) 



with 



C = V (det 7 - det a) - d (det 7 - det a) - det Y 



We have numerically evaluated the quantity 



Lfb = max L . 



(37) 



(38) 



by using Eas. (|3l7l) - (|331) . In doing so we have accounted 
for the stability condition M > and for the Heisenberg 
uncertainty condition T + (i/2)íl > [l6|. Here Í7 is the 
Standard symplectic form 



íï = 



üj 
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1 

-1 



(39) 



As matter of fact, not all real vàlues of A are admitted 
since the addition of feedback may lead to unstable or 
unphysical states. Thus, the maximization (|38|l is per- 
formed on a 'vàlid' range of A's vàlues. The extension of 
such a range clearly enlarges as r\ decreases. Any way, the 
maximum is always achieved for vàlues —(1/2) < A < 0. 

The quantity Lfb of Eq.JJHJl is shown in Fig|3as func- 
tion of x for diffcrent vàlues of 77. Notice that Lfb for 
ï] = (lower curve) corresponds to the case of no feed- 
back. It telis us that in absence of feedback, entangle- 
ment arises for x > and icreases with it till the value 
L = 1. Instead, much higher vàlues can be obtained 
by the feedback action, even with low vàlues of rj. The 
benefit of feedback increases by increasing rj and it is 
particularly manifest for high vàlues of x- I n the limit 
case of -q — * 1 by approaching the instability x ~ * 1/2 
entanglement increases indefinitely. This is because in 
such a case the feedback is able to completely recycle the 
information lost by the system into environment, i.e. it 
somehow suppresses the amplitude damping, making a 
maximally entangled state achievable [2J, [lÇ| . 

It is worth noting that the presented scheme can be 
implemented by using an experimental set up similar to 
that of Ref.|l7|, where the two fields emerging from an 
optical parametric oscillator were separately subjected to 
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ing of a suitable quadrature combination [20j , which then 
might play the role of X\ — X2 ■ 

In conclusion, we have shown the possibility to greatly 
improve the steady state entanglement in an open quan- 
tum system by using a feedback action. In the stud- 
ied case, where Gaussian states were involved, Marko- 
vian (direct) feedback sufhces to reach the goal, while we 
guess that other feedback procedures, like state estima- 
tion based feedback |2l| , could be more powerful in other 
contexts. 



FIG. 2: Logarithmic negativity L/t as function of X- Curves 
from bottom to top are for 77 = 0, 0.3, 0.5, 0.7, 0.99. 



homodyne measurements and the relative currents com- 
bined. Then one has to employ driving field modula- 
tors to realize the feedback action. Any delays in clas- 
sical communication (including the feedback loop) must 
be much smaller than the typical time scale of the sys- 
tem and much smaller than the inverse of relevant 
bandwidth. Experimcntalists can certainly achievc loop 
delays smaller than 10 -8 s 0] which already satisfies the 
above requirements for set ups similar to |l7J • From this 
point of view Markovian (direct) feedback has also the 
advantage of not requiring information processing aftcr 
measurement. To give a numerical example, in Ref.|l7l| 
an overall efficiency of ~ 0.7 was achieved, which already 
means the possibility of ~ 250% improvement of entan- 
glement close to the threshold (see Fig|3Jl. 

The feedback action results a Gaussian operation from 
the quadratic form of Ea. ((T7fl together with Eci. ((T3|l . 
Moreover, the presence of vacuum noises Wj in such 
equations legitimate us to intend the feedback action as 
local operation supplcmentcd by classical communication 
(LOCC). As matter of fact, in absence of interaction 
(x = 0) the feedback action is unable to generate en- 
tanglement as it is evident from FigEl However, in Refs. 
[ï^| it was shown the impossibility to enhance (distill) en- 
tanglement by means of Gaussian LOCC. The key point 
is that, in contrast with Refs. 0], we have considered 
continuous measurement and continuous communication 
between the two parties, including the exchange of noise 
of quantum origin. And the correlations among Wj and 
Xj, yj are able to reduce the quantum noisy effects on 
the system, hence to enforce the interaction between the 
subsystems. This corresponds to an overall nonsepara- 
ble Gaussian map in the sense of Ref. 0], though the 
real physical operations are LOCC. Thus, the presented 
approach may shed further light on the subject of entan- 
glement distillation. 

It is also to remark that the our model extends behind 
Eq.(HJ to any other quadratic interaction, because Gaus- 
sian state's entanglement is always related to the squeez- 
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